February 7, 2008 7:37 Proceedings Trim Size: 9in x 6in 



noon2004 



o 
o 

CN ■ MINIBOONE STATUS 

3 



< 



X 



X 



J. L. RAAF* FOR THE BOONE COLLABORATION 

University of Cincinnati 
Department of Physics, M.L. 0011 
Cincinnati, OH 45221, USA 
E-mail: jlraaf@fnal.gov 



> 

oo 

o 

O . 

00 , MiniBooNE, the Mini Booster Neutrino Experiment at Fermilab, will confirm or 

C _~i i refute the existence of the neutrino oscillation signal seen by the Liquid Scintillator 

Neutrino Detector (LSND) Experiment at Los Alamos National Laboratory. The 
C ~ ~) ' experiment will search for the appearance of electron neutrinos in a beam of muon 

neutrinos. This work presents preliminary results from the first round of analyses of 
charged current quasi-elastic events, neutral current n° events, and neutral current 
elastic events. The neutrino oscillation analysis is not presented in this work; it 
O l' is a blind analysis which will not be presented until the full set of data has been 

, collected. 
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1. Introduction 



'_ ' MiniBooNE, the Mini Booster Neutrino Experiment at Fermilab, has been 



running smoothly and taking data for over a year now. The experiment 
intends to make a definitive statement about the neutrino oscillation signal 
seen by the Liquid Scintillator Neutrino Detector (LSND) Experiment at 
Los Alamos National Laboratory 1 . Since both solar and atmospheric neu- 
trino oscillations have been well-established in recent years, confirmation of 
the existence of another oscillation signal with its own distinct m 2 splitting 
would indicate a need for further extensions to the Standard Model in order 
to accommodate this signal. 

The Fermilab Booster operates on a 15 Hz cycle, sending between 2 
and 3 Hz of its total to MiniBooNE. This 8 GeV proton beam from the 
Booster is directed into a 71 cm long beryllium target located in a magnetic 
focusing horn. The horn sends positively charged particles into a 50 meter 
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decay region. Any charged particles that do not decay in the 50 meters are 
absorbed by or decay in the absorber located at the downstream end of the 
decay region. The neutrinos produced in the charged particle decays then 
travel through approximately 500 meters of dirt to reach the detector. By 
the time the beam reaches the detector it is very pure - almost entirely made 
up of muon neutrinos. Only about 0.5% of the beam is electron neutrinos. 
The average energy of a neutrino entering the MiniBooNE detector is 800 
MeV. 

The detector consists of a sphere within a sphere, where the inner signal 
region is optically isolated from the outer veto shell 2 . The entire detector is 
instrumented with 8-inch Hamamatsu photomultiplier tubes (PMTs): 1280 
in the main signal region of the detector (10% coverage) and 240 in the veto 
region 3 . The entire detector is filled with approximately 800 tons of pure 
mineral oil 4 . 

Beam arrives at the target in 1.6 [is bursts. The data acquisition system 
(DAQ) receives a signal from the Booster indicating that beam will be 
arriving imminently. This signal triggers the DAQ to open a window and 
record all detector activity for 19.2 fis, starting 4.8 ^s before beam arrives. 
Fig. 1 shows the distribution of events during the beam window. The left 
panel shows the 1.6 (is wide peak of the beam particles arriving in the 
middle of the beam window. 




Average Tank Time ( sec) Average Tank Time ( u sec) Average Tank Time ( u sec) 

Figure 1. The left panel shows recorded tank activity during 19.2 /js beam window 
with no cuts. The center panel adds a simple veto cut. The right panel adds a further 
requirement of greater than 200 hits in the main volume of the detector. 



With the application of only two very simple cuts, the background on 
the distribution can be drastically reduced. Requiring no more than 6 hits 
in the veto region removes events with a coincident cosmic muon occurring 
in the beam window, shown in the center panel of Fig. 1. The further 
requirement of at least 200 hits in the main tank region, shown in the 
right panel, largely removes the chance of an electron from muon decay 
occurring during the beam window. With only these two cuts, the signal 
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to background ratio is greater than 1000:1. As further confirmation that 
the particles observed are actually coming from the beam, we check the 
reconstructed direction of these particles with respect to the beam direction. 
The majority of the particles align with the beam direction. 



2. Neutrino Interaction Analyses 

The three main analyses that are currently underway at MiniBooNE are 
charged current quasi-elastic scattering (CCQE), neutral current tt° pro- 
duction (NC 7T°), and neutral current elastic scattering (NCE). At Mini- 
BooNE energies, about 39% of neutrino interactions are due to CCQE, 7% 
NC 7T°, and 17% NCE. From the dataset of about 1 x 10 20 protons on target 
analyzed, we already have an abundance of events to study for each of these 
analyses. 



2.1. Charged Current Quasi- Elastic Events 
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Figure 2. Reconstructed neutrino energy, E„ , (left) and negative square of the four- 
momentum transfer, Q 2 , between the neutrino and the muon (right) for charged 
current quasi-elastic events. The black dots indicate data, the grey rectangles indicate 
the Monte Carlo predictions including the current level of systematic uncertainties. All 
distributions arc normalized to unit area. 



Since charged current quasi-elastic interactions are among the most 
prevalent events at MiniBooNE energies, there are already more than 30,000 
events in the study sample. It is important to understand these interactions 
because v e CCQE events with similar kinematics and cross section will be 
the main signal for the — ► v e oscillation search, and these interactions 
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will help to constrain the intrinsic v e flux in the beam. In addition, com- 
parison of the actual CCQE cross section with the expected cross section 
will allow MiniBooNE to make a measurement of disappearance in this 
channel for Am 2 <~ 0.1 — 10 cV 2 . The CCQE event selection requires a series 
of cuts to obtain a high-purity sample of well-reconstructed v„ events. The 
cosmic ray-induced activity is reduced by the method discussed in section 
1. Events in which the neutrino energy cannot be reliably reconstructed are 
rejected via a fiducial volume cut. Finally, the selection requires a single 
muon-like Cercnkov ring and scintillation light production consistent with 
a v^n —> fi~p interaction. 

The CCQE analysis 5 determines the neutrino energy, E Ul using the 
outgoing muon energy, E^, and outgoing muon angle relative to the neu- 
trino beam direction, cos 9^. Reconstructed neutrino energy for data is 
compared to Monte Carlo predictions, both normalized to unit area, in the 
left panel of Fig. 2. The error bars on this figure include our current as- 
sessment of the effect of the major sources of systematic uncertainties. In 
particular, systematic errors associated with the flux, CCQE cross 
section, and the properties of light production and transmission in the Mini- 
BooNE detector are taken into account. Substantial reductions in these 
systematics are expected as the analysis progresses. 

The right panel of Fig. 2 shows the measured and expected distribu- 
tions in Q 2 , the negative of the squared four-momentum transfer between 
the incoming neutrino and the outgoing muon in the two-body CCQE scat- 
tering process. The behavior of the Q 2 distribution at low Q 2 values is 
particularly important in experimentally validating how well nuclear ef- 
fects are modeled in the Monte Carlo description of neutrino interactions 
on the bound target nucleons of the MiniBooNE detector. 

2.2. Neutral Current tt° Events 

Understanding the rate and kinematics of ir° production in neutrino neutral 
current interactions is critical for MiniBooNE, since these events are a major 
background for the — > v e appearance search. The preliminary analysis 
of the mass spectrum and kinematics of 7r°'s produced in neutral current 
interactions observed in MiniBooNE is presented here. 

Fig. 3 shows the reconstructed invariant mass of beam triggers which 
satisfy the NC ir° event selection criteria. As with the CCQE analysis, 
cosmic ray-induced activity during the beam trigger is eliminated by the 
method discussed in section 1. Further cuts on fiducial volume and en- 
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Mass = 0.1363 ± 0.0008 GeV/c: 
Width = 0.0192 ± 0.0008 GeV/[T 
Nuln. lt° = 2509 ± 1 13 Events . 
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Figure 3. Reconstructed invariant mass of beam triggers (open circles with statistical 
error bars). Fitted shapes are Monte Carlo-based parameterizations of the contribution 
from background (dashed curve) and signal (NC resonant and coherent 7r°'s, solid minus 
dashed curve). Note that some of the background events contain 7r°'s, so the peak near 
m^o in the dashed curve is expected. 
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Figure 4. Extracted n° yields in p n o (top), energy asymmetry (middle), and cos 9^o 
(bottom) bins comparing Monte Carlo and data, each normalized to unit area. Error 
bars on data are statistical only. Monte Carlo error bars are estimated systematic errors 
from the yield extraction procedure only, added in quadrature with statistical errors. 
Statistical errors alone are indicated by the innermost hashes on MC error bars. 



ergy of gammas from ir° decay are applied to ensure a good reconstruction 
of the event. The fitted curves are Monte Carlo-based parameterizations, 
where the signal contribution arises from NC resonant and coherent single 
7T° events, and the background contribution comes from all other events. 
The background shape is determined by Monte Carlo simulations. A back- 
ground peak near m^o is expected, since the background does contain some 
7r°'s produced in final state interactions and some 7r°'s from multi-pion 
events. For approximately 1 x 10 20 protons on target, we extract 2425 ±107 
signal-like NC ir° events in the fit to data in the preliminary analysis. The 
fitted mass peak agrees well with the nominal ir mass. 
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The number of neutral current 7r°'s seen in data is also extracted in 
bins of calculated variables of interest (p^o , cos O^a , ^qr|N ) to produce a 
distribution of each variable with the binned yields. The signal fraction in 
each bin is extracted via a fit to the invariant mass plot for events in that 
bin. The distributions for each variable are compared with expectations 
from Monte Carlo simulations, both distributions normalized to unit area. 

Fig. 4 (left) shows the tt° momentum distribution for data and Monte 
Carlo. For all parts of this figure, the error bars shown on data are sta- 
tistical only. Monte Carlo errors are statistical added in quadrature with 
estimated systematic errors stemming from the yield extraction procedure. 
The fraction of 7r°'s reconstructed at high momentum drops off mainly due 
to the fall-off in the neutrino flux spectrum. 

The energy asymmetry of the 7's from w° decay is shown in Fig. 4 (mid- 
dle). It falls off just below 0.9 because of the minimum energy requirement 
for each of the 7's in the decay; events where one or both of the 7's have less 
than 40 MeV are cut. The angular distribution of 7T°'S relative to the beam 
direction is shown in Fig. 4 (right). The distribution is sensitive to the pro- 
duction mechanism - coherent tt°s are much more strongly forward-peaked 
than resonant 7r°'s. This distribution may help us to determine how much 
coherent production contributes to the overall rate of NC tt° production. 

2.3. Neutral Current Elastic Scattering Events 

MiniBooNE can also easily isolate neutral current clastic scattering 
events 6 with only a few simple cuts. The primary background to NCE 
scattering events is cosmic muons, and the subsequent Michel electrons 
from the muon decay. Random strobe triggers, which contain the beam- 
unrelated backgrounds from cosmics and environmental radioactivity, are 
used to make a "beam-off" subtraction. The method is equivalent to that 
described in section 1, but the cuts are optimized for very low energy events. 
This is shown in Fig. 5. In the left panel, note the 1.6-/is-wide peak due to 
beam neutrino interactions. The falling exponential is caused by the decay 
of cosmic muons. Once the subtraction is performed, a flat distribution, 
except for the expected beam spike, remains. 

Neutral current elastic interactions should produce relatively low energy 
nucleons, and therefore a low number of tank hits. They are also expected 
to leave a relatively high amount of scintillation light in the detector, since 
only nuclear fragments are detected. Fig. 6 shows the strobe-subtracted 
distributions of tank hits for events with scintillation fraction > 0.5. The 
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Figure 5. Average event time strobe subtraction. In the left panel, the 1.6 /ts wide 
peak located at 4.6 - 6.2 /is is from beam neutrinos. The center panel shows random 
strobe triggers. The right panel shows the beam window with random strobe triggers 
subtracted off. 




Figure 6. Tank hits with strobe subtraction compared to relatively normalized Monte 
Carlo for events with scintillation fraction > 0.5. 



Monte Carlo is relatively normalized to the data for shape comparison. 
The histogram shows the relatively high fraction of neutral current events 
expected in this sample. 



3. Expectations for the Future 

MiniBooNE has been running smoothly for more than a year now. As 
of February, 2004, the experiment has collected over 200,000 contained 
neutrino candidates with the 1.8 x 10 20 protons on target received. We 
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expect to have 5 x 10 20 protons on target by June, 2005. The three non- 
oscillation analyses described in this paper, CCQE, NC 7r°, and NCE, are 
well underway. The main items remaining to be studied and understood for 
these are the systematic errors associated with the flux, cross sections, 
and the properties of light production and transmission in the MiniBooNE 
detector. These analyses are necessary for determining the two oscillation 
results, disappearance and v e appearance, which are the main goals of 
this experiment. 
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